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Abstract. Using a sample of 407 RR Lyrae stars extracted from SDSS Stripe 82, we 
study the degree of substructure in the Galactic halo. We identify overdensities associated 
with the known substructures of the Hercules-Aquila Cloud and the Sagittarius Stream, 
and find a further previously-unknown substructure, the Pisces Overdensity, at ~ lOOkpc 
from the Sun. Together, the three substructures account for ~ 80% of our RR Lyrae sam- 
ple. We also study the density distribution of RR Lyraes in the halo and find that the profile 
is best fit by a broken power law with an inner slope of -2.4 and a break radius of ~ 25 kpc, 
consistent with previous studies. We conclude that the halo is predominantly made up of 
the debris from in-falling satellites, with little or no underlying smooth component. 

1 Introduction 

Simulations of the formation of our Galaxy suggest that its halo must have formed, either entirely (e.g. 
O)) or in part (e.g. |2|), via a series of minor mergers; small subhalos caught by the potential of the 
Milky Way (MW) fall in and are tidally stripped as they do so, leaving streams of stars along their 
orbits. These streams and the remnants of the progenitor satellites might be seen as overdensities in 
the halo, telling the story of the formation history of the Galaxy. Here we study substructure using RR 
Lyrae stars in a ~ 300-square-degree patch of sky known as Stripe 82. 



RR Lyraes are variable stars that are commonly used as tracers of substructure (e.g. |I31|4]|5|); their 
distinctive lightcurves make them easy to identify and their well-defined absolute magnitudes allow for 
reliable distance estimates. Furthermore, they are believed to be a common phase of stellar evolution, 
so they are ubiquitous, and they are bright, so that they can be detected out to large distances. 

Stripe 82 is an equatorial strip of the SDSS footprint in the Southern hemipshere, where the observ- 
ing strategy was to observe a small area many times. The result is a data set (e.g. |6 1) with an eight-year 
baseline for which we have an average of 30 epochs per object (and as many as 80 epochs in the most 
extreme cases), making it well suited for identifying and studying variable stars. Furthermore, there 
is data in five photometric bands, important for classifying and extracting stellar populations, and a 
subset of the data have spectra and, hence, metallicity estimates. 

Bramich et al. (|7|) realised the huge value of such a data set and created two catalogues in order to 
best exploit the data. They are the Light-Motion-Curve Catalogue (LMCC), which contains lightcurves 
for all objects in the stripe, and the Higher-Level Catalogue (HLC), which contains over 200 derived 
quantities, e.g. mean magnitudes, mean colours and variability statistics. 

Stripe 82 is of further interest because the Hercules-Aquila Cloud (|8|) is known to extend into the 
area covered by the stripe and the orbit of the Sagittarius Stream (e.g. |9 |) also crosses the stripe. 

We selected a sample of 407 RR Lyraes in Stripe 82 using the LMCC and HLC via a series of 
colour, period and metallicity cuts. We then used the information in the catalogues to calculate absolute 
magnitudes for our sample and then distances accurate to ~ 5%. Further details on the selection method 
can be found in 1 10|. Here, we present the main results from the analysis of our RR Lyrae sample. 
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Fig. 1. The distribution of RR Lyraes in SDSS Stripe 82 in iieliocentric distance and Riglit Ascension. Tlie 
distribution is not uniform in Right Ascension and does not fall off smoothly with distance, suggesting that there 
is a lot of substructure in this part of the halo. Three distinct substructures can be identified: the Hercules-Aquila 
Cloud (highlighted in red) and the Sagittarius Stream (highlighted in orange) were both previously known to exist. 
The Pisces Overdensity (highlighted in blue) was identified during the course of this study. 



2 Substructure in the halo 

The final sample of RR Lyraes is shown as function of heliocentric distance and Right Ascension 
(RA) in Figure [T] along with the distance errors. The distribution is not smooth in RA, neither does 
the number of RR Lyraes fall off smoothly with distance; clearly the halo is far from uniform in this 
part of the sky. 

We are able to identify three significant substructures in our distribution (for details of member 
selection see |10|). There is a large overdensity at high RA and small distance, highlighted in red, 
that has both distance and RA consistent with the Hercules-Aquila Cloud. This substructure alone 
contains 237 stars, nearly 60% of our sample. In this range of distance and RA, we would expect to 
find a smooth bulge component present in our sample as well as RR Lyraes belonging to the halo. 
We attempted to separate the two samples using metallicities, however they proved impossible to 
disentangle. We found that our metal-rich, intermediate metallicity and metal-poor samples all had the 
same spatial distribution (see [iOJ, Figure 15). 

At low RA and small distance, we find a second clump of RR Lyraes which have positions and 
distances consistent with the Sagittarius Stream where it crosses the stripe. Our Sagittarius sample, 
highlighted in orange, contains 55 stars, which is a further 15% of the total sample. 

The third overdensity - the Pisces Overdensity - is highlighted in blue. It is more distant at ~100 
kpc and contains 28 stars, or 7% of our sample; it looks as if the structure could be bifurcated, although 
with so few stars, it is not possible to be certain. This structure was also identified by Sesar et al. (PTTIl. 
"Structure J") although they did not confirm whether it was real or simply an artefact in their data. 

First we ask: is it real? The error bars at these distances are large enough to be visible in the figure 
but not large enough to explain the presence of this substructure. We also analysed extinction changes 
across the stripe but found very little variation and conclude that Pisces is indeed a real structure. 

Models of the Sagittarius Stream do not predict that the stream would be detected in this part of 
the sky, so it is unlikely that this clump is associated. The Magellanic Stream (MS) has been detected 
in this part of the sky, but is found at smaller distances. Furthermore, the MS has only been traced in 
gas and is very diffuse in this region; if we were to detect the first stars in the MS then it would be 
more likely in a denser region. So we are unable to associate Pisces with any known substructure. 

We can make some simple comparisons with other known objects in the halo in order to learn 
more about the Pisces Overdensity; although, of course, any comparisons that we make can only tell 
us about the part of Pisces that we have detected in the stripe. First we compare Pisces with MW 
dwarf spheroidal (dSph) galaxy Carina ( [,12,1 ). which also hes at a heliocentric distance of ~100 kpc. 
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Fig. 2. Left: the fraction of accessible volume probed by the SDSS in Stripe 82 as a function of Galactocentric 
distance. It is possible to detect RR Lyraes out to 100 kpc and the brightest RR Lyraes out to 130 kpc. Right: the 
cumulative number of RR Lyraes with Galactocentric distance; this distribution is best fit by a broken power law 
(red line). We also see a sharp drop off beyond ~40 kpc, though it is not clear whether this is a real property of 
the halo or due to the large amount of substructure present in Stripe 82. 



Assuming similar stellar populations and detection efficiencies and using simple scaling arguments, 
we estimate a mass of ~ lO^M© for Pisces in the stripe. 

We can also compare with our Hercules- Aquila sample, which has the advantage of being detected 
via the same method, but this cloud is less well understood than Carina. Again, we assume similar 
stellar populations to calculate a mass of ~ 10*Mq for Pisces in the stripe, reassuringly similar to our 
previous estimate. We also find that Pisces is 4-5 times more diffuse than Hercules-Aquila in the stripe. 

Since we completed this work, there has been follow-up spectroscopy for a handful of stars in 
the Pisces Overdensity. KoUmeier et al. {(13]) obtained spectra for eight stars, of which five had very 
similar velocities. The dispersion of these stars suggests that Pisces is a disrupted dSph galaxy and not 
a disrupted globular cluster (GC), although whether the remnant is bound or unbound is unclear. The 
remaining three stars from their study were clearly offset from the main population but had similar 
velocities to each other, which suggests that there may indeed be a bifurcation, or that we could be 
looking at two different objects. Sesar et al. ([14]) obtained spectra for five stars, of which one was 
also in the Kollmeier et al. (|13|) sample, bringing the total of Pisces stars with spectra to 12. They 
also found evidence for two kinematic groups and concluded that Pisces is spatially extended and is 
likely to be unbound. 



3 Density distribution of RR Lyraes in the halo 

In order to study the number density of RR Lyraes in the halo, we first need to consider the fraction 
of accessible volume probed by the SDSS in Stripe 82 as a function of Galactocentric distance, which 
we show in the left panel of Figure|2] We generated RR Lyrae populations, assuming a Gaussian lumi- 
nosity function with a mean and dispersion the same as that found for our RR Lyraes. The magnitude 
limits were set to be the same as those defining our sample. We see that the survey reaches out to at 
least 100 kpc before it starts to decline, and the brightest (Mz ~ 0.1) RR Lyraes are visible all the way 
out to 130 kpc. 

In the right panel of Figure |2] we show the cumulative number of RR Lyraes with Galactocentric 
distance, where we have adjusted by the fraction of Galactocentric volume sampled and assumed a 
selection efficiency of 1. Smooth halo models often invoke density power laws with indices of -3, 
but we find the best fit to the data is a double power law, which we show in red. Such models were 
first proposed by Saha (|15|), who noticed that the density fell off more rapidly beyond 25 kpc; our 
estimate of the turn-over radius of 23 kpc is in excellent agreement. And our estimate of -2.4 of the 
inner slope is also consistent with a study by Miceli et al. (| 16 |) of a very large sample of RR Lyraes 
out to 30 kpc. 
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We see a sharp drop off in the number density of RR Lyraes beyond ~40 kpc. The idea of such 
a drop off in RR Lyrae number is certainly not new, but it has not be wholly agreed upon. Ivezic et 
al. (|3|) found evidence for an "edge" in a study of RR Lyraes found in SDSS commissioning data; 
although, they later applied the same method to the a larger sample that extended out to 70 kpc and 
found no such edge anymore. Vivas & Zinn (|5|) also conducted a similar study and found no edge 
out to 60 kpc. Our data suggests that "edge" is perhaps too strong a term for the phenomenon, as the 
number of RR Lyraes declines sharply rather than being completely cut off, but there is definitely a 
drop off and it is real. It cannot be due to the properties of the survey - remember, we are able to see 
even the faintest RR Lyraes out to 100 kpc. 

Of course, we saw in the previous section that there is a lot of substructure in the stripe, so its 
arguable how much we can really say about the density of the halo from this data. The Hercules- 
Aquila Cloud alone accounts for almost 60% of our RR Lyraes, and all 3 substructures contain nearly 
80% of the RR Lyraes. Hercules-Aquila and Sagittarius dominate the stripe out to ~40 kpc, so it could 
be their extents that explain the double-power law that we find. Possibly there is little or no underlying 
smooth component in the halo and we are seeing mostly substructure. 



4 Summary 

We extracted a set of RR Lyrae stars from the SDSS DR6 data of the region known as Stripe 82 and 
used them to study substructure in the halo. We recovered two substructures that had been previously 
identified, and found a third that had been suggested, though not confirmed. This Pisces Overdensity 
is found at ~100 kpc and cannot be connected to any known substructures. The Magellanic Stream 
is also found in that part of the sky but it has only been traced in gas and is found at much smaller 
distances. The size of the errors and extinction changes across the stripe cannot explain it away. 

We also looked at the distribution of RR Lyraes in the halo. We found that RR Lyraes in Stripe 
82 are visible out to 100 kpc, and that the brightest are visible all the way out to 130 kpc. We also 
found that the density distribution is best fit by a broken power law, with a sharp fall off at ~40 kpc, 
consistent with previous results. However, we should bear in mind the dominance of substructure in the 
stripe that might naturally explain these results. Certainly our RR Lyrae sample suggests that, if there 
is a smooth background halo component underneath the complexity left behind by satellite accretions, 
then it is only minimal. 

Acknowledgements: I am grateful to Wyn Evans and Vasily Belokurov for their invaluable assistance 
with this project. I also acknowledge funding from the Science and Technology Facihties Council. 



References 

1. Bullock, J.S. and Johnston, K.V., ApJ 635, (2005) 931-949 

2. Abadi, M.G., NavaiTO, J.F. and Steinmetz, M., MNRAS 365, (2006) 747-758 

3. Ivezic, Z. et al., AJ 120, (2000) 963-977 

4. Vivas, A.K. et al., AJ 127, (2004) 1 158-1175 

5. Vivas, A.K. and Zinn, R., AJ 132, (2006) 714-728 

6. Adelman-McCarthy, J.K. et al., ApJS 175, (2008) 297-313 

7. Bramich, D.M. et al., MNRAS 386, (2008) 887-902 

8. Belokurov, V. et al., ApJL 657, (2007) L89-L92 

9. Belokurov, V. et al., ApJL 642, (2006) L137-L140 

10. Watkins, L.L. et al., MNRAS 398, (2009) 1757-1770 

11. Sesar, B. et al., AJ 134, (2007) 2236-2251 

12. DairOra, M. et al., AJ 126, (2003) 197-217 

13. Kollmeier, J.A. et al., ApJL 705, (2009) L158-L162 

14. Sesar, B. et al., ApJ 717, (2010) 133-139 

15. Saha, A., ApJ 289, (1985)310-319 

16. Miceli, A. et al., ApJ 678, (2008) 865-887 



